load("data/RData/all.RData")
load("data/RData/hama.RData")
library(tidyverse)
library(summarytools)
library(mapview)上位25%以上の密集地帯・ホットスポットです。
map6<-map5 %>%
mutate(hotspot=if_else(exp2>=5.8,1,0)) %>%
filter(hotspot==1)
sfa<-intersec %>%
select(geometry,total) %>%
mutate(hot=if_else(total>=27,1,0)) %>%
filter(hot==1)
mp<-mapview::mapview(sfa,
zcol="hot",
col.regions = c("orange"),
map.types =c("Esri.WorldImagery","OpenStreetMap"))
mapview::mapview(x=map6,
zcol="hotspot",
col.regions = "red",
color="transparent",
map.types =c("OpenStreetMap","Esri.WorldImagery"))+
(mp)mapview::mapview(x=map5,
zcol="exp2",
color="transparent")